Infrared cutoff dependence of the critical flavor number in three-dimensional QED 
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We solve, analytically and numerically, a gap equation in parity invariant QED3 in the presence 
of an infrared cutoff fj, and derive an expression for the critical fermion number N c as a function of 
fx. We argue that this dependence of N c on the infrared scale might solve the discrepancy between 
continuum Schwinger-Dyson equations studies and lattice simulations of QED3. 
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Parity invariant quantum electrodynamics in 2+1 di- 
mensions with N flavors of four-component massless 
fermions (QED 3 ) has been attracting a lot of interest 
for almost two decades. While it was often regarded as a 
nice polygon for studying nonperturbative phenomena in 
gauge field theories, such as dynamical mass generation, 
recently the model has found applications in condensed 
matter physics, in particular in high-T c superconductiv- 
ity 0. 

QED3 is ultraviolet finite and has a built-in intrinsic 
mass scale given by the dimcnsionful gauge coupling e, or 
a = e 2 N/8, which plays a role similar to the A scale pa- 
rameter in QCD. In the leading order in 1/N expansion, 
it was found that at large momenta (p 3> a) the effec- 
tive coupling between fermions and gauge bosons van- 
ishes (asymptotic freedom) whereas it has a finite value 
or infrared stable fixed point at p <C a 3] . This behavior 
was shown to be robust against the introduction of higher 
order 1/N corrections (J ||. Studies of the gap equa- 
tion for a fermion dynamical mass in the leading order 
in the 1/N expansion have shown that massless QED3 
exhibits chiral symmetry breaking whenever the num- 
ber of fermion species N is less than some critical value 
N c , which is estimated to be in the region 3 < N c < 5 
(N c = 32 /ir 2 ~ 3.2 in the simplest ladder approximation 
la). A renormalization group analysis gives 3 < N c < 4 
7J. Below such a critical N c , the U(2N) flavor sym- 
metry is broken down to U(N) x U(N), the fermions 
acquire a dynamical mass, and 2A^ 2 Goldstone bosons 
appear; for N > N c the particle spectrum of the model 
consists of interacting massless fermions and a photon. 
In addition, it was argued that the dynamical symme- 
try breaking phase transition at N — N c is a conformal 
phase transition 0, . The last one is characterized by 
a scaling function for the dynamical fermion mass with 
an essential singularity 0, @ • 

On the other hand, there is an argument due to Ap- 
pclquist et al. [lOj that N c < 3/2. The argument is based 
on the inequality /ir < /tjv where / is the thermody- 
namic free energy which is estimated in both infrared 
and ultraviolet regimes by counting massless degrees of 
freedom. 

The above-described version of QED 3 , the so-called 
parity invariant noncompact version, was studied on a 



lattice 11]. Recent lattice simulations of QED3 with 
N > 2 have found no decisive signal for chiral symme- 
try breaking I n particular, for N = 2, Ref. 01 
reports an upper bound for the chiral condensate (ipip) 
to be of order 5 x 10~ 5 in units of e 4 . We recall that 
for quenched, N = 0, QED 3 both numerical and an- 
alytical 0] studies have shown that chiral symmetry is 
always broken and numerically {ijji/j) ~ 5 x 1CP 3 . The 
above result of Ref. 0] seemingly contradicts studies 
based on Schwinger-Dyson (SD) equations which advo- 
cate dynamical symmetry breaking for N — 2 and seems 
to favor the estimate of Appelquist et al. (we mention 
also that work in progress is being done on lattice sim- 
ulations of QED3 with a single, N — 1, fermion flavor 
0). 

One of the major problems in studying dynamical sym- 
metry breaking using lattice simulations is to obtain con- 
trol of finite size effects, which play a nontrivial role due 
to the presence of a massless photon I n terms of 

the intrinsic lattice spacing a, three dimcnsionless length 
scales appear in lattice simulations: the size of the lattice 
L ~ 10 — 50, the dimcnsionless lattice coupling constant 
(3 oc l/(e 2 a), and the bare fermion mass mga. In order 
to establish dynamical symmetry breaking the role of all 
these scales in the problem should be well under con- 
trol. At present, lattice discretization appears to be well 
understood; however, the finite size effects appear non- 
trivial and are most likely the source of the discrepancy 
between lattice and continuum studies of QED3 |l2| . 

In this paper, we present a possible explanation for the 
discrepancy between recent lattice and continuum studies 
of QED3 by studying analytically and numerically the 
gap equation in the presence of the infrared (IR) cutoff 
fx, where /i is inversely related to the size L of a lattice. 
Since the characteristic ratio e 2 //z for continuum QED3 
turns out to be proportional to the ratio L/0 in lattice 
simulations 0], a comparison between two approaches 
can be made. We will show that the presence of an IR 
cutoff reduces the value of the critical number N c and 
derive the relationship between N c and fj,. Recently, a 
similar gap equation, but with a massive photon, was 
studied numerically in Refs. |l6l IT^ in connection with 
applications of QED3 to high-T c superconductors (for an 
earlier numerical study in QED3, see Ref. |l8j). 
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The gap equation of massless QED3 with IR cutoff /1 
reads (compare with Eq. (3.3) of Ref.Q]) 



dk- 



fcE(fc) k+p + a 
! + S 2 (fc) n |jfc-p| +a' ' 1 



where A = 8/Nir 2 and the Landau gauge is used. Equa- 
tion is the simplest approximation to the SD equa- 
tion for the fermion self-energy which neglects corrections 
to the fermion wave-function renormalization and vertex 
corrections. Extensive studies showed that the nature 
of chiral symmetry breaking which emerges from Eq. 
with fi = remains qualitatively unchanged after includ- 
ing higher order corrections. In what follows we solve 
Eq. numerically but in order to be able to treat it 
analytically one needs to make further approximations. 
Because the integrand is damped at large momenta of in- 
tegration, the main contribution comes from the region 
with momenta k <C a; thus, expanding the logarithm we 
come to the simplified gap equation 



S(p) 



dk 



fcS(fc) 
k 2 + S 2 (fc) 



mhx(k,p). 



(2) 



The scale \x can be identified with the inverse size of the 
lattice in the temporal direction as /1 ~ ir/La since 
both scales fi and it /La represent a minimal fermion mo- 
mentum in continuum and lattice theories, respectively 
\p is related to the photon mass m a (fj, — m 2 /a) in the 
model of Ref. 0]]- Although this argument for the iden- 
tification of [i and it / La (up to a factor of order 1) is 
certainly heuristic, it is not unreasonable. 

On the other hand, the ultraviolet (UV) cutoff a in 
Eq. J2J is m f ac t a physical scale which separates a non- 
perturbative dynamics at k <C a from a perturbative one, 
with momenta much larger than a, where the dimen- 
sionless running coupling a(k) — a/(k + a) [j| is weak. 
Neglecting those perturbative contributions in Eq. 
results in an effective UV cutoff a in Eq. J3J. 

We shall solve the above equation J2J using the 
well-established linearized approximation when the 
momentum-dependent dynamical mass in the denomina- 
tor k 2 +T, 2 (k) is replaced by the constant dynamical mass 
at the lower limit: S 2 (fc) — > E 2 ,, So = E(//). This approx- 
imation is known to work well, especially near the phase 
transition point where So <C ol. As we show, the lin- 
earized equation with an IR cutoff can be dealt with 
analytically and an exact expression for N c as a function 
of a/fj, can be derived, following the method developed 
in Refs. [H HI 

In the linearized approximation, Eq. J5J can be reduced 
to a hypergeometric differential equation 

«(1 - u)S"(«) + ~(1 - u)E'(u) - ~E(u) = 0, (3) 
with the IR and UV boundary conditions 



,dS 
dp 



= 0, 



p— + S 

dp 



= 0, 



(4) 



where u = — p 2 /Sp. The general solution of Eq. © can 
be written as 



S(P) 
So 

where we choose 



ClMl(p) + c 2 u 2 (p), 



,'1 + iv \-iv 3 p 2 

u iW = ^ ( — ;— . — r-; o5 -^2 j . 



(5) 



(6) 



"2(P) 

x F ( - 



p_ 

So 

1 — if \-\-iv 



-{l+iv)/2 



W Sn \ 

1 + Y ;--|)+C.C., (7) 



as a particular pair of independent solutions of Eq. J2J), 
with v — V4A — 1, and where F is the hypergeometric 
function. If the infrared scale is set equal to zero (/i = 0) , 
then the IR boundary condition gives the constraint c 2 = 
0, since the function u 2 is too singular at p = 0. Then, the 
UV boundary condition fixes the relationship between 
the dynamical mass So, f, and a. For nonzero /x, the 
boundary conditions Q determine the mass spectrum, 
giving rise to the equation 



/(So/a, n/a, v) = A X B 2 - A 2 B 1 = 0, 



where 
At 



dui 
dp 



B, 



p—a 



2 dUi 

^ dp 



(8) 



(9) 



p=p 



By making use of the formulas for differentiating the hy- 
pergeometric function p3| the real functions At and Sj 
are expressed as 



Ai = F 
A 2 -- 



1 + iv 1 — iv 1 a 2 



(10) 



So\ 



1 — if 1 + if 
x F 1 — ,—— :1 



+ c.c. 



(11) 



and 
Si 



/i 2 (1 + f 2 ) 



S 2 12 



x F 1 



1 + if 1 — if 5 /i 2 
A + ; 2 ;_ SI 



(12) 



\ -(l+»i/)/2 /, . • x 

*--"<£) ^ 



x F 1 



1 + if 1 — if 



1 + 

* 2 ' M 2 

x [Re(F(f))cos/?2 + Im(F(f))sin/3 2 ], 



(13) 
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where /3 2 = {v/2) ln(/i/So) — tan 1 v 1 and F(y) is short- 
hand notation for the hypergeometric function: 

( 1 + if 1 — if if Sn\ 

F ^ = F { 1+ —--—-> 1+ r>-^)- ^ 

Since we study the critical behavior, we can always as- 
sume that So <SC a, and therefore we can use the asymp- 
totic expressions for A\ and A 2 : 



Ax « ^ 



a 
So 



and 



where 



a 

s7 



.4, 



-1/2 



-1/2 



r(l/2)r(ii//2) / a 



iv/2 



ra((l + ii/)/4) VSo 
|c(f)|\/l + i^ 2 cos (a 2 + 



c.c. 



/ \ -1/2 



COS «2, 



(15) 



(16) 



a 2 = - In =r + tan 1 v, = argc(z/), 



c( V ) 



r(3/2)r(w/2) 



r((l + ii/)/4)r((5 + ii/)/4) 



By making use of (2.10.2) of [23|, we can write Bi in a 
form similar to -B 2 : 

I U \ (1 + If) , 



Re(F(i/)) 



x cos(/3 2 + 6») + Im(F(z/)) sin(/3 2 + i 



(17) 



Finally, using the above expressions (fT^|l - fT7| . the gap 
equation JHJ can be expressed as 



Re(F(z/)) 



x sin(/3 2 - a 2 ) - Im(F(z/)) cos(/3 2 - a 2 ) 



(18) 



One can convince oneself that for /i < So the last equa- 
tion is reduced to 



f a 
cos — In — 

2 So 



tan 1 v I = 0, 



(19) 



which gives the well-known result for the dynamical mass 
exhibiting the essential singularity at v — > 0: 



Sq = aexp 



2 (7m + tt/2 -9- tan" 1 v) 



(20) 



We recall that only the solution with n = 1 corresponds 
to the stable ground state. Note also that for N close to 
iV c (v ~ 0) the dynamically generated mass So is much 
less than the scale a, providing a hierarchy of mass scales 
in the model under consideration. 

On the other hand, for So <C /i, we can use the power 
expansion of F(v) to get the equation for the dynamical 
mass S . By expanding F{y) of Eq. for S <C fi, we 
obtain 

S 2 

Re(F(i>)) = 1 + -I- p(f 2 ) cos <j>, (21) 
M 



Im(F(z/)) = -§p(i/ 2 ) sin ( 



where 



1 {2h + v 2 ){l + v 2 ) 



tan 



-i ^(13 



4 + v 2 



10 - Iv 1 



Thus Eq. ifTBl) reduces to 



M 2 



sin ( £ In ^ 



2 tan 



1 

/ , (^ 2 )sin(§ln^ + 2tan 



-l 



f 



v - 



(22) 

(23) 
(24) 

(25) 



Since, for < v < 1 0, is negative, we find that a non- 
trivial or real solution for So arises when v exceeds the 
critical value v c determined by the equation 



v c , a. 
— In — 

2 /! 



2 tan 1 v c = 7r, 



(26) 



so that the sin in the numerator of Eq. I|25(l becomes 
negative (for v > v c ). Note that the form of the gap 
equation l|18|l is different in two regions (i <C So and 
/i 3> So: while in the first one (/i <C So) we observe 
oscillations in the mass variable, in the second one (/i ^> 
So) such oscillations disappear. This is reflected also in 
the character of the mass dependence on the coupling 
constant [compare Eqs. (|20|l and l|29ll below]. In general, 
it can be shown that Eq. i|18|) has n nontrivial solutions 
where the number n is given by 



v a 
— In — 

27T /! 



tan v 



(27) 



and the symbol [C] denotes the integer part of the num- 
ber C. For small f c , Eq. I|26|l gives 

(28) 



2 + (1/2) ln(a//i) ' 

A similar functional dependence of N c was guessed in 
Ref. by fitting their numerical study of the gap equa- 
tion with a massive photon. 

Near the critical point, we find from Eq. (|25(l the mean- 
field scaling law for So, 



M 2 



a{v - v c ), 



(29) 
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FIG. 1: The numerical solution of Eq. (0 and the analytical 
solution of Eq. @ given by Eq. Il2(i> for N c versus the ratio 
ire 2 /fi or the lattice ratio L//3. 



where the positive "nonuniversal" constant of propor- 
tionality, <x, is 



1 



(it — 2 tan 1 v c ) 



1 



1 



sin(7r + (f> c ) ' 



with <p c the value of at v = v c . 

We now compare the analytically obtained critical 
value of v c (Eq. (|26pl or iV c as a function of the ratio 
7re 2 /n (a = Ne 2 /8) with the numerical solution of the 
integral equation Q which was obtained by using nu- 
merical methods described in Ref. [24]] ■ Basically, start- 
ing with some initial guess for S(p), the integral equa- 
tion is iterated using splines on a logarithmic momentum 
scale, until sufficient precision is reached. Both the an- 
alytical solution for Eq. (J2J) and numerical solution for 
Eq. are depicted in Fig. ^ The analytical solution 
of Eq. agrees reasonable well with the numerical so- 
lution of Eq. QJ for values of N > 1. However, Eq. (J2J 
(hence Eq. provides a rather poor description of 

dynamical symmetry breaking for small values of N ( 



actually it does not allow quenched limit in contrast to 
Eq.QJ). The bending of the curve for the analytic solu- 
tion for < N < 1 is due to the linear dependence of the 
effective UV cutoff a on N. Furthermore, the numerical 
solution shows the existence of a critical ratio 7re 2 / ^ w 40 
for N = (i.e., quenched QED3) below which no chiral 
symmetry breaking occurs. 

From Fig. ^ it can be extracted that in order to find 
chiral symmetry breaking for N = 2 at least a ratio 
we 2 /n w 5 x 10 3 is required. Note that in order to draw 
decisive conclusions on chiral symmetry breaking on a 
lattice, the volume of the lattice L 3 must be large not 
just in dimensionless units but also compared to any dy- 
namically generated correlations in the system. In partic- 
ular, the physical volume (£//3) 3 is required to be large 
fl2| . Typical lattice simulations use at the most a ra- 
tio L//3 w 50 — 100 and appear to lie outside the region 
for dynamical symmetry breaking. Indeed, by associat- 
ing La ~ tt//j, and e 2 ~ l/(/3a), we have L/(3 ~ ire //z; 
thus, dynamical symmetry breaking near N — 2 requires 
a ratio L/f3 w 5 x 10 3 . This explains the absence of any 
signs of chiral symmetry breaking in lattice simulations 
with N — 2. For N = 1 the ratio is L//3 w 240 for the 
numerical solution of Eq. QJ; therefore, we believe that 
lattice simulations with a single fermion flavor are crucial 
in establishing whether there is a critical number N c in 
QED3 and how this value depends on the size L of the 
lattice. Moreover, as is evident from Eq. I|29|l . the scaling 
near the critical point is expected to be of the mean-held 
type and such a scaling can also be checked using lattice 
simulations. 

We are grateful to V.A. Miransky for fruitful discus- 
sions and comments and S. Hands for valuable corre- 
spondence. We thank A.H. Hams for help in the nu- 
merical programs. V.P.G. acknowledges support from 
the Natural Sciences and Engineering Research Council 
of Canada. The research of V.P.G. has been also sup- 
ported in part by the SCOPES-projects 7 IP 062607 and 
7 UKPJ062150.00/1 of Swiss NSF and by NSF Grant No. 
PHY-0070986. 



* On leave of absence from Bogolyubov Institute for The- 
oretical Physics, 03143, Kiev, Ukraine. 

[1] R. D. Pisarski, Phys. Rev. D 29, 2423 (1984). 

[2] M. Franz and Z. Tesanovic, Phys. Rev. Lett. 87, 257003 
(2001); Z. Tesanovic, O. Vafek, and M. Franz, Phys. 
Rev. B 65, 180511 (2002); M. Franz, Z. Tesanovic, and 
O. Vafek, ibid. 66, 054535 (2002); I. Herbut, ibid. 66, 
094504 (2002); ibid. Phys. Rev. B 66, 094512 (2002). 

[3] T. Appelquist, M. Bowick, D. Karabali, and L. C. R. Wi- 
jewardhana, Phys. Rev. D 33, 3704 (1986). 

[4] D. Nash, Phys. Rev. Lett. 62, 3024 (1989). 

[5] V. P. Gusynin, A. H. Hams, and M. Reenders, Phys. 
Rev. D 63, 045025 (2001). 

[6] T. Appelquist, D. Nash, and L. C. R. Wijewardhana, 



Phys. Rev. Lett. 60, 2575 (1988). 
[7] K.-I. Kubota and H. Terao, Prog. Theor. Phys. 105, 809 
(2001). 

[8] V. A. Miransky and K. Yamawaki, Phys. Rev. D 55, 
5051 (1997). 

[9] V. P. Gusynin, V. A. Miransky, and A. V. Shpagin, Phys. 

Rev. D 58, 085023 (1998). 
[10] T. W. Appelquist, A. G. Cohen, and M. Schmaltz, Phys. 

Rev. D 60, 045003 (1999). 
[11] E. Dagotto, A. Kocic, and J. B. Kogut, Phys. Rev. Lett. 

62, 1083 (1989); Nucl. Phys. B334, 279 (1990); V. Az- 

coiti and X.-Q. Luo, Mod. Phys. Lett. A 8, 3635 (1993). 
[12] S. J. Hands, J. B. Kogut, and C. G. Strouthos, Nucl. 

Phys. B645, 321 (2002), S. J. Hands, J. B. Kogut, 



5 



L. Scorzato, and C. G. Strouthos, Nucl. Phys. B (Proc. 
Suppl.) 119, 974 (2003). 
[13] S. J. Hands and J. B. Kogut, Nucl. Phys. B335, 455 
(1990). 

[14] K. Kondo and H. Nakatani, Mod. Phys. Lett. A 5, 407 
(1990); V. P. Gusynin, A. W. Schreiber, T. Sizer, and 
A. G. Williams, Phys. Rev. D 6 0, 065007 (1999). 

[15] P. Maris and D. Lee, hep-lat/0209052 D. Lee and 
P. Maris, Phys. Rev. D 67, 076002 (2003). 

[16] G. Liu and G. Cheng, |con£ma t/0208061 Phys. Rev. D 
67, 065010 (2003). 

[17] T. Pereg-Barnea and M. Franz, Phys. Rev. B 67, 
060503(R) (2003). 

[18] K. Kondo and H. Nakatani, Prog. Theor. Phys. 87, 193 
(1992). 



[19] Because of the antiperiodic boundary condition for 
fermions in the temporal direction, the coefficient is n 
instead of 2n. 

[20] V. A. Miransky, V. P. Gusynin, and Yu. A. Sitenko, Phys. 

Lett. 100B, 157 (1981). 
[21] P. I. Fomin, V. P. Gusynin, V. A. Miransky, and 

Yu. A. Sitenko, Riv. Nuovo Cimento 6(5), 1 (1983). 
[22] V. Gusynin, M. Hashimoto, M. Tanabashi, and K. Ya- 

mawaki, Phys. Rev. D 65, 116008 (2002). 
[23] H. Bateman and A. Erdelyi, Higher Transcendental 

Functions (McGraw-Hill, New York, 1953), Vol. 1. 
[24] V. P. Gusynin, A. H. Hams, and M. Reenders, Phys. 

Rev. D 53, 2227 (1996). 



